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Abstract 

We use configuration space methods to write down one-dimensional integral rep- 
resentations for one- and two-loop sunrise diagrams (also called Bessel moments) 
which we use to numerically check on the correctness of the second order differential 
equations for one- and two-loop sunrise diagrams that have recently been discussed 
in the literature. 
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1 Introduction 



Sunrise-type diagrams have been under investigation since many years. Exact analytical 
results can be obtained only for special mass or kinematic configurations such as for the 
equal or zero mass cases or for the threshold region. For example, threshold expansions 
of the non-degenerate massive two-loop sunrise diagram have been studied in Refs. [TJ [2]. 
The construction of differential equations for the corresponding correlator function provides 
some hope that by solving these differential equations, a general analytical solution can be 
obtained. Recently, mathematical methods were used to construct the coefficients of such a 
differential equation in a systematic way [3] . This work supplements the work of Kotikov [I] 
and Remiddi et al. jSJ El [7] on the same subject. While traditionally the correlator is 
calculated in momentum space, configuration space techniques allow for a surprisingly 
simple solution for sunrise-type diagrams: The correlator in configuration space is just a 
product of single propagators which in turn can be expressed by modified Bessel functions 
of the second kind. Transforming back to momentum space, one ends up with a one- 
dimensional integral over Bessel functions, known as Bessel moments [SJ [H] . As outlined in 
a series of papers [101 Ell EE2l 131 [2j [HI [151 [16], the corresponding one-dimensional integral 
can be easily integrated numerically for an arbitrary number of propagators with different 
masses in any space-time dimension. Therefore, configuration space techniques can be 
used to numerically check the differential equations for the correlator function obtained by 
other means. This will be detailed in this note. 

The paper is organized as follows: In Sec. 2 we introduce the configuration space 
techniques which will be used in Sec. 3 to check the differential equations for one-loop 
sunrise-type diagrams. In Sec. 4 we check the differential equations for the two-loop sunrise 
diagrams for the equal mass case, while in Sec. 5 we will deal with nondegenerate masses. 
Our conclusions can be found in Sec. 6. Even though the configuration space techniques 
are well suited to treat general space-time dimensions, we will mainly deal with the 

case of D = 2 space-time dimensions in this paper. For reasons of simplicity, throughout 
this paper we work in the Euclidean domain. The transition to the Minkowskian domain 
can be obtained as usual by a Wick rotation (or, equivalently, by replacing p 2 — > — p 2 ). 
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2 Configuration space techniques 



n( x m ) - L / c " * - fVH^ /ox 



In configuration space, the n-loop n-particle irreducible correlation function 

n(x) = (0|Tj(x)j(0)|0) (1) 
connecting the space-time points and x is given by the product of the propagators, 

n+1 

IL(x) = HD(x, mi ), (2) 

where the free propagator of a particle with mass m in D-dimensional (Euclidean) space- 
time is given by 

1 r e t(j}X ^d D p ( y mx) x K\(mx) 
(2^p J p 2 + m 2 = (2tt) a+1 x 2 

(D = 2 A + 2). K\(z) is the McDonald function (modified Bessel function of the second 

kind). Note that p and x in the integral expression of Eq. ([3]) are D-dimensional Lorentz 

vectors, i.e. p ■ x = p^x^ - , while the quantity x in the rightmost expression of Eq. ([3]) (and, 

therefore, also in the argument of the propagator) denotes the absolute value x = ^fx^x^ . 

In the limit mx — > at fixed x, the propagator simplifies to 

D M ) = ^f^ = i (4) 
l ' J (2vr)W p 2 An^x 2 ^ l J 

where T(A) is Euler's Gamma function. 

In this note we write the n-particle irreducible correlator function in (Euclidean) mo- 
mentum space. The momentum space n-particle irreducible correlator function is given by 
the Fourier transform of the n-particle irreducible correlator function U(x) in configuration 
space, 

f[(p) = / U(x)e i{px) d D x. (5) 

As a product of propagators, U(x) in (j5J) depends only on the absolute value x = ^x^x^, 
Therefore, one proceeds by first integrating the exponential factor over the D — 1 dimen- 
sional hypersphere. We write d D x = x D ~ 1 d D xdx where dPx denotes the D — 1 dimensional 
integration measure over the D — 1 dimensional hypersphere. The integration of the expo- 
nential factor over the unit sphere gives 

e«r-*)d D x = 2n x+1 (^~* 



(p*)d D x = 2n x+1 (^) " J x (px). (6) 
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J\(z) is the Bessel function of the first kind. As before, p and x on the right hand side of 
Eq. (jSJ) stand for the absolute values p = yjp^p 1 and x = tJx^x^. Therefore, the correlator 
in momentum space depends only on the absolute value of the momentum, 



This is the central formula for our numerical verification of the correctness of the differential 
equations. 

3 The one-loop case 

In Ref. [5], Remiddi explains how to obtain the differential equation for the one- loop sunrise- 
type diagram with arbitrary masses and dimensions. By applying the integration-by-parts 
technique to the correlator in momentum space, recurrence relations can be obtained. 
Finally, Euler's theorem for homogeneous functions connects the loose ends of the iterative 
steps involving partial derivative with respect to p 2 . We have numerically checked all these 
steps and have found numerical consistency - up to Stokes' contributions due to surface 
terms in integer space-times dimensions. 
To be more precise, the integral 



for v = k,p (or a linear combination of both) leads to a surface term which can be assumed 
to vanish (up to Stokes' contributions). The integration- by-parts technique consists in 
calculating the integral explicitly and one then expresses the result in terms of scalar 
integrals 
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The (two) resulting recurrence relations read 












+ a 2 {{p 2 + ml + ml)S(ati, a 2 + 1) — S(a 1 — 1, a 2 + 1) — S(a>i, a 2 )j , (10) 

-a% [{p 2 -ml + ml)S(ax + 1, a 2 ) — S(ax + 1, a 2 — 1) + S(ax, a 2 )J 

+ a 2 {{p 2 + m\ - ml)S(ax, a 2 + l) — S(ax — 1, a 2 + 1) + S(ax, a 2 )) • (11) 
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Eq. (Ilip can be replaced by Eq. (fTOI) with the two lines interchanged, 

= D S(ai, a 2 ) + 2a 2 (mlS(ai,a 2 + 1) — S(cti, a 2 ) 

+ «i [{p 2 + m\ + m|)5'(ai + 1, a 2 ) - 5(«i + 1, a 2 - 1) - 5(cm, a 2 )) • (12) 

It is obvious that Eq. ( ITT]) is reproduced as difference of Eq. (TTOl) and Eq. (fT2|) . Therefore, 
one has to check only Eq. (TTUj) . We will perform this numerical check for the parameter 
choice a.\ — a 2 — 1 which is relevant for the differential equation, and for D = 2 space- 
time dimensions. As mentioned in Ref. [3J, even though other dimensions are feasible, this 
choice avoids singular contributions and serves for the simplest integrand. The equation 
to be checked is 

2m 2 5(2, 1) + (p 2 + m\ + m|.)5(l, 2) = 5(1, 1) + 5(0, 2). (13) 

Starting from 

/d 2 k 1 1 r°° 

j2^ 2 (k 2 + rn 2 )((p-k) 2 + rn 2 ) = 2^Jo ^^oK^oM^, (14) 

the integrals S(a\, ot 2 ) with higher (integer) values of can be obtained as partial deriva- 
tives with respect to the masses, 

2midmi ' dm 2 ' dm 2 J (2n) 2 (k 2 + m 2 ) ((p — k) 2 + m 2 ) 

d k 1 

: 5(2,1) (15) 



and accordingly 



In addition one has 



(2n) 2 (k 2 + mj) 2 ((p - k) 2 + m 2 2 ) 



' ° 5(1,1) = 5(1,2). (16) 



5(2,0) 



2m 2 dm 2 
d 2 k 



(2tt) 2 (fc 2 + ml) 2 47rmf 



The derivative can be expressed by K' Q (z) = —Ki(z). The higher order integrals in the 

configuration space representation are given by 

-1 d 1 
--—5 1,1 = — 2 

lvti\ om\ Airmf Jo 



5(2,1) = — — a _ 5(1, 1) = — — 2 / Jo(px)(mix)Ki(mix)K (m 2 x)x dx, 



-Id 1 r°° 

5(1,2) = - — - — 5(1,1) = K I Jo(px)K (m 1 x)(m 2 x)Ki(m 2 x)xdx. (It 

2m 2 am 2 Anmi Jo 



5 



Therefore, Eq. ( TT3|) in the configuration space representation reads 



J (pa;) 



2(m2x) 2 (mix)Ki(mix)K (m2x) 



+ ((P x ) 2 + (mix) 2 + (m2x) 2 )K (mix)(m2x)Ki(rri2x) 

dx 



2(rri2x) K (mix)Ko(rri2x) 



x 



(19) 



We were able to check this equation numerically for different values of p as function of 
mi and m 2 . The 3D-plot in MATHEMATICA shows stochastic fluctuations around the 
exprected value of 1 of the order of 10~ 9 . 

Euler's theorem of homogeneous functions leads to the differential equation 



Because of J'q(z) = —Ji(z), Eq. (120]) can be translated to 



(20) 







{pxy 



(px)Ji(px)KQ(mix)K (m2x) + Jo(px)(mix)Ki(mix)K (m2x) 



+ J (px)K (mix)(m 2 x)Ki(m2x) - 2J (px)K (m 1 x)K (m2x) 



— (21) 

x 



We checked on the latter relation with an even better precision of the order of 10 -13 . 
The differential equation 

d 

(p 2 + (mi + m 2 ) 2 ) (p 2 + (mi - m 2 ) 2 ) q^S(1, 1) = -{p 2 + m\ + m 2 2 )S{l, 1) 



(22) 



in Ref. j5] is obtained by inserting the recurrence relations into Euler's differential equa- 
tion (|20|) . Using the configuration space representation, Eq. fT22|) reads 



1 



— ((px) 2 + (mix + m 2 x) 2 ^ ({px) 2 + {m\X — m 2 x) 2 ^) (px) J±(px) 



2(px) 

+ ((P x ) 2 + (rriix) 2 + (m 2 x) 2 ) J (px) 



K (mix) K (m 2 x) — . 

x 



(23) 



This equation could be checked with a precision of the order of 10 -8 

While Euler's differential equation can be derived from general principles also for the 
configuration space representation, the recurrence relations can be derived only via the 
momentum space representation. If one does not use this technique, it remains unclear 
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why such integral identities exist for general parameters p, m\ and 777.2. In order to check 
whether one can derive further relations by using integral identities in configuration space, 
we have used the general ansatz 



Jo(px) 



A (px, mix, m2x)Ko(mix)Ko( y m2x) 



+ Ai(px, mix, m2x)(mix)Ki(mix)Ko(m2x) 
+ A 2 (px, mix, m2x)K (mix)(m2x)K 1 (m2x) 



dx 

x 



(24) 



where m 1; m 2 ) = a i0 p 2 + a iX m\ + aj 2 m| (i = 0, 1, 2). By chosing random values for p, 
mi and m 2 , and solving the resulting system of equations, one obtains 

a o = 0, a i = —2a, a 2 = — 2 + 2a, 



a^Q — a, c^ii — O; O12 — 2 — a, 



^20 = 1 — a, a 2 i = 1 + a, a 2 2 = 1 — a, 



(25) 



where a is an arbitrary parameter. This leads to Eq. (T191) and 



JoGra) 



2 ((m 2 x) 2 — (mix) 2 ) K Q {m l x)K Q {m2x) 
+ ((P^) 2 + {rnix) 2 — (m 2 x) 2 ) (mix) K^rriix) K (rri2x) 
— (ipx) 2 — (mix) 2 + (m 2 x) 2 ) Ko(mix)(m2x)Ki(m2x) 



dx 

x 



(26) 



which is the difference of Eq. (Tl9|) and the same equation with mi and m 2 interchanged. 
We conclude that no more recurrence relations can be found that go beyond Eq. ( Fl9l) . 



4 The two-loop case with equal masses 

The differential equation for the two-loop degenerate sunrise diagram has been given in 
Ref. [7j. It reads 

(2p 2 (p 2 + m 2 )(p 2 + 9m 2 ) (4^\ + (3(4 - D)p A + 10(6 - D)m 2 p 2 + 9£>m 4 ) 



+ (D - 3) ((£> - 4)p 2 - (£> + 4)m 2 ) 5(1, 1, 1 



(D — 4) 2 7T 2 ' 



(27) 



which simplifies to 



W + *W + ^) f - V + (*« + 20mV+9m 4 )-^ + (p* + W)Wl, 1,1) 3 



<ip 2 



<ip 2 



(28) 

in £) = 2 space-time dimensions. We write 5(1, 1, 1) in a form which is easily adapted to 
the non-degenerate mass case to be discussed later on. One has 



5(1,1,1) 



(2tt)- jo 



Jo(px) K (mix) K (m 2 x) K (m 3 x)x dx . 



(29) 



Differentiation of Eq. ( 1291) gives 

dp 2 ' ' 2pdp 

(px) Ji(px) K (mix) Ko(rri2x) Ko(rri3x)x dx 



-1 



-5(1,1,1) 



2(2?r)V Jo 



(30) 



and 



(px)J 1 (px) + -(px) 2 (J 2 (px) - J {px)) 



4(2vr) V Jo 
x Ko(mix)K (m2x)K (m3x)x dx 



(31) 



Returning to the degenerate mass case the differential equation (|28p can be translated to 

1 



3 
2 



(p 2 + m 2 )(p 2 + 9m 2 ) ( (px)Ji(px) + -(px) 2 (^(px) — Jo(px)) 



4p 2 Jo 

- 2(3p 4 + 20m 2 p 2 + 9m 4 )(px)J 1 (px) + Ap 2 (p 2 + 3m 2 )J (px) 



K (mx) 3 xdx. (32) 



where J' Q ' (z) = —J[(z) = {Ji{z) — Jq(z))/2 is used. Because the result contains the second 
derivative of the Bessel function, one can use Bessel's differential equation 



z 2 j';( z ) + z j' x ( z ) + ( z 2 -x 2 )j x ( z ) = o 



(33) 



for A = to compactify the result, 

^2 J (4p 2 (p 2 + 3m 2 ) - (p 2 + m 2 )(p 2 + 9m 2 )(px) 2 ) J (px) 



3 
2 



2p (p + 5m )J\{px) 



Ko(mx) 3 x dx. 



(34) 



These results have been checked with a precision of the order of 10" 



5 The two-loop case with arbitrary masses 



The final section of this paper is devoted to the second order differential equation, derived 
for the two- loop sunrise diagram with arbitrary masses in Ref. [3]. After adjusting the 
normalization, the differential equation can be written as 

2 



po(—p 2 ) ( ~r?) +pi(-p 2 )^ + P2(-p 2 ) 

\ ap A I ap A 



where the coefficients Pi(t) (i = 0, 1, 2, 3) are given by 
Po(t) 



t(t — (mi + m 2 + m 3 ) 2 ) (t - (—mi + m 2 + 

X (t - (mi + m 2 - m 3 ) 2 ) (3t 2 - 2M W0 t - M 200 + 2Mi 



5(1 ' M) = lMp (35) 



m 3 ) 2 ) (t - (mi - m 2 + m 3 ) 2 ) x 
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(36) 



Pi{t) 



9f 



' - 32M 100 t 5 + (37M 200 + 70M no )t 4 - (8M 300 + 56M 210 + 144M m )t 3 
(13M 400 - 36 M 310 + 46M 220 - 124M 2n )t 2 



BM 500 + 24M 4 io - 16M 320 - 96M 3U + 144M 221 )t 
(M 600 - 6M 5 io + 15M 420 - 2OM330 + I8M411 - 12M 32 i - 6M ; 



222y 



(37) 



p 2 (t) = 3t 5 - 7Mi 00 t 4 + (2M 200 + 16M no )t 3 + (6M 300 - UM, 



r 21 )t 2 



- (5M 400 - 8M 310 + 6M 220 - 8M 2n )t 
+ (M 500 - 3M 4 i + 2M 320 + 8M311 - 10M; 



(3J 



Ps(t) 



-18t 4 + 24Mi 00 t 3 + (4M 200 - 40M no )t 2 + (-8M 300 + 8M 210 + 48M m )t 

+ (-2M 400 + 8M 3 io - 12M 220 - 8M211) + 2c(t, m x , m 2 , m 3 ) ln(m 2 //i 2 ) 

+ 2c(t, m 2 , m 3 , mi) ln(m 2 //i 2 ) + 2c(t, m 3 , mi, m 2 ) ln(mg/ /i 2 ), (39) 



and where 



i 2 WAiW r (A 2)f 2 W A 3 ) 



M Ai a 2 a 3 — X!( m i) v' ;t 2; v^; 
are monomial symmetric polynomials in m 2 , m 2 and m| and where 

c(t, mi, m 2 , m 3 ) = f— 2m 2 + m 2 . + m 2 ) t 3 

+ (§m\ — 3m^ — 3m 4 — lm\m\ — lm\m\ + Mm 2 .?/! 2 ) t 2 
+ ( - 6mi + 3m^ + 3mjj + \\m\m\ + llm 4 m 3 



(40) 
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+ (2m\ — m\ — m\ — hm\m\ — hm\m\ + m\m\ + m\m\ + Am\m\ + 4m 



+ 3m 4 m 2 + 3m 4 m 4 — §m\m\ + 2m\m\m\ — m\m\m\ — m\m\m^) (41) 



(for details, cf. Ref. [3]). Using Eqs. Q22J, flUD) and flH), one obtains 



p 4 



(px)J 1 (px) + ^(px) 2 (J 2 (px) - J (px)) 



J)l ( — J 

- 2 ^ — (px)Ji(px) + Ap 2 (-p 2 )Jo(px) 



ip 2 (-p 2 ) - ^p-(px) 2 )Jo( P x) 



K (rriix) K (rri2x) Ko(rri3x)x dx 



P 



pz 



Ko(mix) Ko(rri2x) Ko(rn,3x)x dx 



(42) 



where p*(t) = Pi(f) +po(t)/t. Using different values for p and m^, in terms of m-i and 
m 2 we obtain a 3D-plot with MATHEMATICA which shows again stochastic fluctuations 
of the order of 10 -4 . In the course of our numerical checks we were able to identify two 
typos in the coefficients of c(t, mi, m.2, m^) in the preprint version of Ref. [3] which we have 
corrected. 



6 Conclusions 

Using configuration space techniques, we were able to check numerically the differential 
equations for sunrise-type diagrams found in the literature. The precision of our numerical 
test is still quite moderate, but gives sufficient confidence in the validity of the differen- 
tial equations derived by other means. For example, the introduction of artificial "typos" 
in the coefficients of the differential equations are easily discovered. More rigorous tests 
would require the use of more stable integration routines than those provided by MATHE- 
MATICA. For the future we hope to find independent routes to discover further relations 
between Bessel moments which may lead to generalizations of the present findings to cases 
involving three-loop or even higher order sunrise-type diagrams. 
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